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In this proceeding, results of searches for CP violation in charm decays using the full BABAR. dataset are discussed. 
The parameter Aqp in the decay — > K^n^ 1 is determined to be (—0.39 ± 0.13 ± 0.10)%. Meaurements of 
CP violation using T-odd correlations in the four-body decays D+ — > A'+ K^tt + tt~ ' and Df K+K°n+n- 
O ■ " * 10 " 3 ^ ( ^ ± 3.4 (8yst) ) x 10-3, respectively . 

(N 

^ 1. Introduction 

m ' 

In the Standard Model (SM), CP violation (CPV) arises from the complex phase of the CKM quark-mixing 
matrix [![. Measurements of the CPV asymmetries in the K and B meson systems are consistent with expec- 
r— i tations based on the SM and, together with theoretical inputs, lead to the determination of the parameters of 
the CKM matrix. CPV has not yet been observed in the charm sector, where the theoretical predictions based 
^ . on the SM for CPV asymmetries are at the level of 1CP 3 or below Q. An observation of CP asymmetries at 
the level of one percent or greater would be a clear indication of new physics. 



2. Search for CP Violation in the decay D+ K° s iz + |0l 



BABAR searched for CPV in the decay D^ 1 — > K^ir^ by measuring the parameter Acp defined as: 
Si a _ r(D + ^K°n+)-T(D-^K s n-) 

CP r(D+ -^k°tt+) + t(d- ->.fc:°7r-)' [) 

where T is the partial decay width for this decay. This decay mode has been chosen because of its clean 
experimental signature. Although direct CP violation due to interference between Cabibbo-allowed and doubly 
Cabibbo-suppressed amplitudes is predicted to be negligible within the SM K° — K° mixing induces a 
time-integrated CP violating asymmetry of (—0.332 ± 0.006) % Q. Contributions from non-SM processes may 
reduce the value of the measured Aqp or enhance it up to the level of one percent [3, Q • Therefore, a significant 
deviation of the Acp measurement from pure K° — K° mixing effects would be evidence for the presence of 
new physics beyond the SM. Due to the smallness of the expected value, this measurement requires a large data 
sample and precise control of the systematic uncertainties. Previous measurements of Aqp have been reported by 
the CLEO-c ((-0.6±1.0(stat)±0.3(syst))% [|) and Belle collaborations ((-0.71±0.19(stat)±0.20(syst))% 0)- 

We select D ± — > Ks n± decays by combining a K® candidate reconstructed in the decay mode — > ir + Tr~ 
with a charged pion candidate. A K® candidate is reconstructed from two oppositely charged tracks with an 
invariant mass within ± 10 MeV/c 2 of the nominal K® mass §. To obtain the final candidate events, a Boosted 
Decision Tree (BDT) algorithm Q is constructed from seven discriminating variables for each candidate: the 
measured proper decay time t(D ± ), the decay distance in the transverse plane L xy (D ± ), the CM momentum 
magnitude p*(£> ± ), the momentum magnitudes and transverse components with respect to the beam axis for 
both the K® and pion candidates. 

A binned maximum likelihood (ML) fit to the m(K gir^) distribution for the retained Z? ± candidates is used to 
extract the signal yield. The total probability distribution function (PDF) is the sum of signal and background 
components. The signal PDF is modeled as a sum of three Gaussian functions, the first two of them with 
common mean. The background PDF is taken as a sum of two components: a background from Djr — > K^K^, 
where the is misidentified as tt^, and a combinatorial background from other sources. The data and the 
fit are shown in Fig. [TJ All of the fit parameters are extracted from the fit to the data sample apart from 
the normalization of the background due to Df — > K^K^, which is fixed to the value predicted by the MC 
simulation. We determine Aqp by measuring the signal yield asymmetry A defined as: 

A= N D D + + +N D D -> (2) 

where N D +(N D -) is the number of fitted D + — > K®ir + (D~ — > K®ir~) decays. The quantity A is the result of 
two other contributions in addition to Aqp- There is a physics component due to the forward-backward (FB) 
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Figure 1: Invariant mass distribution for K^ ± candidates in the data (black points). The solid curve shows the fit to 
the data. The dashed line is the sum of all backgrounds, while the dotted line is combinatorial background only. The 
vertical scale of the plot is logarithmic. 

asymmetry (Ap p) in e + e~ — > cc, arising from j*-Z° interference and high order QED processes in e + e~ — > cc. 
This asymmetry will create a difference in the number of reconstructed D + and D~ decays due to the FB 
detection asymmetries arising from the boost of the center-of-mass (CM) system relative to the laboratory 
frame. There is also a detector-induced component due to the difference in the reconstruction efficiencies of 
D + K ( Jir + and D~ — > K®tt~ generated by differences in the track reconstruction and identification efficiencies 
for 7r + and it~ . While App is measured together with Acp using the selected dataset, we correct the dataset 
itself for the reconstruction and identification effects using control data sets. BABAR developed a data-driven 
method to determine the charge asymmetry in track reconstruction as a function of the magnitude of the track 
momentum and its polar angle which is shown along with the associated errors in Fig. [2j 

Neglecting the second-order terms that contain the product of Aqp and App, the resulting asymmetry can 
be expressed simply as the sum of the two. The parameter Acp is independent of kinematic variables, while 
App is an odd function of cos6* D , where 0* D is the polar angle of the .D candidate momentum in the e + e~ CM 
frame. If we compute A{-\-\ cos0£,|) for the D candidates in a positive cos9* D bin and A(— | cosf^Q for the 
candidates in its negative counterpart, the contribution to the two asymmetries from Acp is the same, while 
the contribution from App has the same magnitude but opposite sign. Therefore Aqp and App can be written 
as a function of | cos#|j| as follows: 

AM\«»P D \) = A (+\^h\)-AHcos6 h \) (3) 

and 

A CP (\ cos^l) = A(+|C0S ^ I) ; A(HC0S ^ I) . (4) 

The selected sample is divided into ten subsamples corresponding to ten cos 9* D bins of equal width and a 
simultaneous binned ML fit is performed on the invariant mass distributions of D + and D~ candidates for each 
subsample to extract the signal yield asymmetries. Using the asymmetry measurements in five positive and in 
five negative cos9* D bins, we obtain five App and five Aqp values. As Acp does not depend upon cosflj,, we 
compute a central value of this parameter using a y 2 minimization to a constant. The Acp and App values 
are shown in Fig. [3J together with the central value and ±1 a confidence interval for Aqp- We determine Aqp 
to be: 

A cp = (-0.39 ± 0.13 ± 0.10)% (5) 



where the first error is statistical and the second systematic. 



Proceedings of the DPF-2011 Conference, Providence, RI, August 8-13, 2011 



3 




4 5 6 

Pion Momentum (GeV/c) 



© 1 

§0.8 
0.6 
0.4 
0.2 


-0.2 
-0.4 
-0.6 
-0.8 




- 


fc 

0.016" 
o 




0.0 14 "3 




Pi 




0.012 




0.01 




0.008 




0.006 




0.004 


1 


0.002 








t) 1 2 3 4 5 6 

Pion Momentum (GeV/c) 

Figure 2: Map of the ratio between detection efficiency for n + and 7r~ (top) plus the corresponding statistical errors 
(bottom). The map is produced using the numbers of ir~ and ir + tracks in the selected control sample. 
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Figure 3: Acp (top) and Afb (bottom) asymmetries for — > K^tt^ candidates as a function of |cos#jj| in the data 
sample. The solid line represents the central value of Acp and the hatched region is the ±1 a interval, both obtained 
from a \ 2 minimization assuming no dependence on | cos 0* D \ . 

3. Search for CP Violation using T-Odd Correlations in Df g) K°K+n+ir + |fTotl 

A search for CP violation in the decays D + — > K + K®ir + Tr~ and Df — > K + K®ir + ir~ using T-odd correlations 
is described here. We define a kinematic triple product that is odd under time reversal using the vector momenta 
of the final state particles in the Dfs rest frame as 



Under the assumption of CPT invariance, T violation is equivalent to CP violation. 
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Figure 4: The K + K®-k + h mass spectrum a) in the D + , and b) in the Z?j" mass region. The curves result from the fits 
described in the text. The distributions of the Pull values are also shown. 



We study the T-odd correlations by measuring the observable expressed in Eq. ([6]) and then evaluating the 
asymmetry 

= T(C T > 0) - T(C T < 0) 
T ~ T(C T > 0) + T{C T < 0) ' 1 ' 

where T is the decay rate for the process under study. The observable denned in Eq. ([7]) can have a non-zero 
value due to final state interactions even if the weak phases are zero [Tlj |. The T-odd asymmetry measured in 
the CP-conjugate decay process, At, is defined as: 

= r(-c r >o)-r(-c r <o) 
T ~ r(-c T >o) + r(-c T <o)' 1 ' 

where Ct = Pk- 1 (Ptt- x Ptt+)- We can then construct: 

•At = 2 ( a t ~ M) , (9) 



which is an asymmetry that characterizes T violation in the weak decay process [12h14| . 

At least four different particles are required in the final state so that the triple product may be defined using 
momentum vectors only [l^. The D meson decays suitable for this analysis method are D + — > K + K2.ir + n~ , 
Df — > K + K®ir + Tr~ and D° — > K + K~ir + ir~ . The search for CP violation using T-odd correlations in D° — > 
K + K~ir + ir~ has recently been carried out by the BABAR Collaboration, and no evidence of CP violation has 
been observed [Ta |. 

The D + and Dj meson decay candidates are reconstructed in the production and decay sequence: 

e + e~ ->• XD+ s) ; T>+ } -> K + K Q s n + it" ; K° -)■ tt+tT , (10) 

using the events with at least five charged particles. To obtain the final set of signal candidates, the p* , the 
difference in vertex probabilities that the parent meson originates from a common vertex and the primary ver- 
tex, and the signed transverse decay length are combined in a likelihood-ratio test. Fig. U shows the resulting 
K + K®ir + Tr~ mass spectra in the D + and Df regions. For each region, the signal is described by the superpo- 
sition of two Gaussian functions with a common mean value. The background is parametrized by a first-order 
polynomial in the D + region, and by a second-order polynomial in the Df region. We extract the integrated 
yields N(D+) = 21210 ± 392 and N(D+) = 29791 ± 337 from the fits, where the uncertainties are statistical 
only. 

We next divide the data sample into four sub-samples depending on Dr s \ charge and whether Ct (Ct) is 
greater or less than zero, and fit the corresponding mass spectra simultaneously to extract the yields and the 
values of the asymmetry parameters At and At- The triple product asymmetries for Cabibbo-suppressed decays 
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D° -> K + K 7T+7T+ (l6j . D+ ->• iT + i<r°7r + 7r and Cabibbo-favored decays D+ -> if+if °7r + 7r~ are summarized 
in Tab. HI The average of the triple product asymmetries is also included in the table 

^T = \{A T + A T ) (11) 

which is not a CP violating parameter but may provide more information on the final-state interactions in these 
decays. 



Table I: Triple-product asymmetries At, At, At, and Et for the Cabibbo-suppressed decays D° — ¥ K + K tt + tt [lfjl ]. 
D + -> K + K°tt + 7v~ [l(J and the Cabibbo-favored decays D+ -> K + K° s -k + -k~ [l(J. The values quoted in units 10" 3 . 



Asymmetry 


D°/D° 


D+/D~ 


D+/DJ 


A T 


-68.5 ± 7.3 ± 5.8 


11.2 ± 14.1 ± 5.7 


-99.2 ± 10.7 ± 8.3 


A t 


-70.5 ± 7.3 ± 3.9 


35.1 ± 14.3 ± 7.2 


-72.1 ± 10.9 ± 10.7 


At 


1.0 ± 5.1 ± 4.4 


-12.0 ± 10.0 ± 4.6 


-13.6 ± 7.7 ± 3.4 


E7 1 


-69.5 ± 6.2 


23.1 ± 11.0 


85.6 ± 10.2 



The final measurements for At in all decays are consistent with zero, however, the values for the T-odd 
asymmetries are considerably larger in D° and decays. The differences in these values for the various 
decays may indicate a difference in the final-state interactions. The final-state interactions may be responsible 



for the hierarchy of lifetimes and branching fractions jl7| . 



4. Conclusion 

Measurements with the final BABAR dataset achieve the precision at the SM prediction for CP violation in 
charm decays. The systematic uncertainties are at the level of the statistical uncertainties. Current and future 
measurements from LHCb, Belle, and SuperB will face the challenge of reducing these systematic uncertainties. 



References 

1 N. Cabibbo, Phys. Rev. Lett. 10, 531 (1963); M. Kobayashi and T. Maskawa, Prog. Theor. Phys. 49, 652 
(1973). 

2 F. Buccella et al, Phys. Rev. D 51, 3478 (1995). 

3 P. del Amo Sanchex et al. [BABAR Collaboration], Phys. Rev. D 83 071103 (2011). 

4 H. J. Lipkin and Z. Xing, Phys. Lett. B 450, 405 (1999). 

5 I. I. Bigi and H. Yamamoto, Phys. Lett. B 349, 363 (1995). 

6 S. Dobbs et al. (CLEO Collaboration), Phys. Rev. D 76, 112001 (2007). 

7 B. R. Ko et al. (Belle collaboration), Phys. Rev. Lett. 104, 181602 (2010). 

8 K. Nakamura et al. (Particle Data Group), J. Phys. G 37, 075021 (2010). 

9 P. Speckmayer, A. Hocker, J. Stelzer and H. Voss, J. Phys. Conf. Ser. 219, 032057 (2010). 

10 P. del Amo Sanchez et al. {BABAR Collaboration), hep-cx/1105.4410 (2011). 

11 I. I. Bigi and H.-B. Li, Int. J. Mod. Phys. A 24S1, 657 (2009). 

12 W. Bensalem, A. Datta, D. London, Phys. Rev. D 66, 094004 (2002). 

13 W. Bensalem, A. Datta, D. London, Phys. Lett. B 538, 309 (2002). 

14 W. Bensalem and D. London, Phys. Rev. D 64, 116003 (2001). 

15 E. Golowich and G. Valencia, Phys. Rev. D 40, 112 (1989). 

16 P. del Amo Sanchez et al. (BABAR Collaboration), Phys. Rev. D 81, 111103 (2010). 

17 M. Gronau and J. Rosner, hep-ph/1107.1232v2 (2011). 



